function f = objfun(c)
    global sigma_l sigma_u
    f = c(1)*sigma_l + c(2)*sigma_u;
end